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We obtain an exact solution for the motion of a particle driven by a spring in a Brownian random- 
force landscape, the Alessandro-Beatrice-Bertotti-Montorsi (ABBM) model. Many experiments on 
quasi-static driving of elastic interfaces (Barkhausen noise in magnets, earthquake statistics, shear 
dynamics of granular matter) exhibit the same universal behavior as this model. It also appears 
as a limit in the field theory of elastic manifolds. Here we discuss predictions of the ABBM model 
for monotonous, but otherwise arbitrary, time-dependent driving. Our main result is an explicit 
formula for the generating functional of particle velocities and positions. We apply this to derive 
the particle-velocity distribution following a quench in the driving velocity. We also obtain the 
joint avalanche size and duration distribution and the mean avalanche shape following a jump in 
the position of the confining spring. Such non- stationary driving is easy to realize in experiments, 
and provides a way to test the ABBM model beyond the stationary, quasi-static regime. We study 
extensions to two elastically coupled layers, and to an elastic interface of internal dimension d, in 
the Brownian force landscape. The effective action of the field theory is equal to the action, up to 
1-loop corrections obtained exactly from a functional determinant. This provides a connection to 
renormalization-group methods. 

PACS numbers: 75.60.Ej, 64.60.av, 05.10.Cc 



I. INTRODUCTION 

The motion of domain walls in soft magnets [H-Q , fluid 
contact lines on a rough surface , or strike-slip faults 
in geophysics @-[9[ can all be described on a mesoscopic 
level as motion of elastic interfaces driven through a dis- 
ordered environment. Their response to external driv- 
ing is not smooth, but exhibits discontinuous jumps or 
avalanches. Physically, these are seen e.g. as pulses of 
Barkhausen noise in magnets [HI El, or slip instabili- 
ties leading to earthquakes on geological faults fHj - llij ). 
While the microscopic details of the dynamics are specific 
to each system, some large-scale features are universal 
[ToT ]. The most prominent example are the exponents of 
the power-law distributions of avalanche sizes (for earth- 
quakes, the well-known Gutenberg-Richter distribution 
[16l - |l8l |) and durations. 

The Alessandro-Beatrice-Bertotti-Montorsi (ABBM) 
model [J is a mean-field model for the dynamics of an 
interface in a disordered medium. It approximates a d- 
dimensional interface in a d + 1-dimensional system, de- 
fined by a height function u(x,t), by a single degree of 
freedom, its averag c height u(t) = ^ / d d x u(x, t). It 
satisfies the equation of motion 

d t u(t) = F (u(t)) - m 2 [u(t) - w(t)} . (1) 

w(t) is the external driving, and F(u) an effective ran- 
dom force, sum of the local pinning forces. In pj, it was 
postulated to be a Gaussian with the correlations of a 
Brownian motion 

[.F(ui) - F(u 2 )] 2 = 2a\ Ul - u 2 \ , (2) 
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where a > characterizes the disorder strength. 

This model has been analyzed in depth for the case of 
a constant driving velocity, i.e. w 

The distribution of avalanche sizes and durations was 
obtained by mapping (fTJ) to a Fokker-Planck equation 0, 
3] . The mean shape of an avalanche was also computed 
using this mapping [22|, [24|, [25| . These results agree well 
with numerous experiments on systems with long-range 
clastic interactions, realized e.g. in certain classes of soft 
magnets, or in geological faults [H, 0, MM M. 

However, long-range-correlated disorder as in ([2]) is a 
priori an unphysical assumption for materials where the 
true microscopic disorder is, by nature, short ranged. 
Hence in realistic systems, it can only arise as a model 
for the effective disorder felt by the interface. This guess, 
originally made by ABBM based on experiments, turns 
out to be very judicious. 

In [2l|, it was shown that the effective disorder for an 
interface with infinite-range elastic interactions is indeed 
given by ©. This led to the wide belief that the ABBM 
model is a universal model for the center-of-mass of an 
interface in dimension d at or above a certain upper crit- 
ical dimension d c depending on the range of the elastic 
interactions in the system [28[ ■ Much of the popularity of 
the ABBM model is owed to this presumed universality. 
However only recently this assumption was proven for 
short-ranged microscopic disorder u sing the Functional 
Renormalization Group (FRG) 1221 . [23j, a method well 
suited to study interfaces (see [29[ for introduction and 
a short review). This proof required quasi-static driv- 
ing w(t) = vt with v = + . Whether this property also 
holds for finite driving velocity v > 0, and in that case up 
to which scale, requires further investigation. The same 
question for non-stationary driving also remains open. 

There are some hints that non-stationary dynam- 
ics may require a different treatment. For example, 
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avalanche size and duration exponents seem to vary over 
the hysteresis loop [30M32"1 | ■ 

Related is the question of static avalanches, i.e. jumps 
in the order parameter of the ground state upon varia- 
tion of an external control parameter, as e.g. the mag- 
netic field. This has been studied for elastic manifolds 
via Functional RG methods [33143511, a nd for spin glasses 
using Replica Symmetry Breaking (3(| H3] • 

In this paper, we discuss the results given by the 
ABBM model when the driving w(t) is a monotonous, 
but otherwise arbitrary function of time. While this 
misses important and interesting physics of AC driving 
and the hysteresis loop [H| , it is much more general than 
the cases treated so far. We will give an analytic solu- 
tion for arbitrary driving, and then specialize to exam- 
ples such as the relaxation of the velocity u(t) after the 
driving is stopped, and the response to finite-size "kicks" 
in the driving force, w(t) = woS(t). This should allow 
to clarify the range of the ABBM universality class by 
comparing these predictions to experiments and further 
theoretical work. Such non-stationary driving can easily 
be realized e.g. in Barkhausen noise experiments, where 
w(t) is the external magnetic field, and can be tuned as 
desired. 

This paper is structured as follows. In section [TT] we 
review the approach to the ABBM model through the 
Martin-Siggia-Rose (MSR) formalism. The MSR formal- 
ism maps disorder averages over solutions of the stochas- 
tic differential equation ([1]) to correlation functions in a 
field theory. In [22|, [23[ this method was used to compute 
the Laplace transform of the p-point probability distri- 
bution of the velocity in the ABBM model, via the so- 
lution of a non linear "instanton" equation. From it, 
the avalanche shape and duration distributions were ob- 
tained for quasi- static driving, in agreement with the re- 
sults of [HI, [25|. Here we extend the method of Ref. 
f22l l23l | and show that it is even more powerful: For any 
monotonous (but not necessarily stationary) driving w(t) 
the resulting field theory can be solved exactly. We give 
an explicit formula for the generating functional of the 
particle velocity it. In scction lTOl wc apply this solution to 
several examples. In particular, we derive the law for the 
decay of the velocity after the driving is stopped, which 
may easily be tested in experiments. In section IIVI we 
extend the method to variants of the ABBM model with 
additional spatial degrees of freedom. This includes the 
generalization of the ABBM model to a (i-dimensional 
interface submitted to a quenched random force with the 
correlations of the Brownian motion, a model whose stat- 
ics was studied in [3~i| . For this more general model, un- 
der monotonous driving, wc show that the action of the 
field theory is not renormalized in any spatial dimension 
d. In section|Vl we compute the generating functional for 
the particle position u, which is more subtle than the one 
for the velocity ii. In sections |VT1 and |VII[ we summarize 
the results and mention possible extensions. In particu- 
lar, wc explain why non-monotonous motion requires a 
separate treatment, and does not follow from the present 



results. 



II. SOLUTION OF THE NON-STATIONARY 
ABBM MODEL 

For understanding the physics of ([I) , one would like to 
know the joint probability distribution for arbitrary sets 
of velocities u(ti)...ii(t n ), averaged over all realizations of 
the random force F. This is encoded in the generating 
functional 

G[A, w] = P*, (3) 

where 7TT denotes disorder averaging. One then recovers 
e.g. the generating function e AM ('°) of the distribution of 
ii(to) by setting A(t) = XSft—to), and similarly for n-time 
correlation functions. 

Our main result is an explicit formula for G in the 
case of monotonous but non-stationary motion. Given 
the distribution of velocities Po(ui) at an initial time ti, 

we claim that G ti := e is 

G ti [X,w] =e m / dViPoMe*™**. (4) 

Jo 

Here u(t) is the solution of an instanton equation (22l.[23|: 

d t u{t) - m 2 u(t) + cru(t) 2 = -X(t). (5) 

Boundary conditions are u(oo) = 0; X(t) is assumed to 
vanish at infinity. Note that u(t) only depends on A(i), 
i.e. the type of observable one is interested in, but not on 
the driving w(t). The latter only enters in 

In the following, we are mostly interested in the case 
when the initial time ti — > —00. Our observables will be 
local in time, so that X(t) decays quickly for t — > ±00. 
Then, u{U) — > and (@| becomes independent of initial 
conditions, 

G[X,w] =e ™ 2 / t ( 6 ) 

To prove ([5]), we first discuss how a closed equation for 
the velocity variable can be formulated. We then use 
the Martin-Siggia-Rose formalism to transform it to a 
field theory, and evaluate the resulting path integral to 
obtain ©. Both steps use crucially the assumption of 
monotonous motion. 



A. Velocity in the ABBM model 

The equation of motion for the velocity u(t) is obtained 
by differentiating ([T]): 

d t u{t) = d t F (u(t)) - m 2 [u(t) - w(t)} . (7) 

A priori, to determine the probability distribution of u(t), 
one needs u(0) and it(0), since the random force depends 
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on the trajectory u(t) and not just on ii. However, un- 
der the assumption that all trajectories are monotonous 
(u(t) > for all times t), the probability distribution of 
ii(t) is independent of u(0). Indeed, under this assump- 
tion, one can replace dtF(u(t)) by a multiplicative Gaus- 
sian noise which only depends on u(t). More precisely, 
we can set d t F(u(t)) - 



vAi(t)C(i) where £(*)£(*) = 
2aS(t — t'). To see this explicitly, consider the generating 
functional 



H[X] 



J t X(t)d t F(u{t)) 

-f tt , \(t)\(t')%d t d t ,\u(t)-u(t')\ 



(8) 



Since u(t) > at all times, we know that [39[ 

d t d v \u{t) - u(t')\ = u(t)d t >sga(u(t) - u(t')) 



u(t)d t >sgn(t - t') = -2u(t)6(t - t'), (9) 



and hence 



H[X] = e J't x( - t)2 ' j{l( - t) = e^Wv™^'. 



(10) 



Note that for monotonous driving, the monotonicity as- 
sumption ii(t) > is enforced automatically if it at t = to 
[H: 



u(t ) > 0, w(t) > for all t > t Q 
u(t) > for all t > t . 



(11) 



In this way we see that for monotonous motion, (J7) is 
a closed stochastic differential equation for the velocity 
ii(t). Given an initial velocity distribution P(u(0)), it 
can be solved without knowledge of the position u(0). 



B. MSR field theory for the ABBM velocity 

The Martin-Siggia-Rosc (MSR) approach allows us to 
express (J3j) , averaged over all realizations of F in j[7| in a 
path integral formalism, following [lE US HI El IS El ■ 

Introducing the Wick- rotated MSR response field u(t) 
and averaging over the disorder, one gets: 

G[X,w] = J Vfaule-^^+St^W, (12) 
S[ii, u]= [ u(t) [d t u{t) + m 2 (ii(t) - w(t))] 
d t d v \u{t)-u{t')\u{t)u{t'). 



Since we consider only paths where ii(t) > at all times, 
using ([9]) we can rewrite the action as 



S[u, u] = j {u(t) [d t u{t) + m 2 (u(t) - w(t))] 

-au{t)u{t) 2 } . (13) 

The key observation which allows to evaluate this exactly 
was first noted in p2l |23| : The action is linear in u(t). 



This means that the path integral over u can be evalu- 
ated, giving a (^-functional. Instead of using this in the 
limit of v — >• as in [22], El , one can write more generally: 

G[X,w] = JV[u,u]e m * /**(*)*(*) x 

xe It u{t)[dtu(t)-m 2 u{t)+ait{t) 2 +\(t)] 

= JV[u] e ™ 2 / t s(t)«>(t) x 

xS [d t u{t) - m 2 u(t) + au(t) 2 + X(t)] , 

This then reduces to ([6]) with u(t) given by ([5]). Note that 
the Jacobian from evaluating the (5-functional is indepen- 
dent of w(t). We assume in the following that for w{t) = 
we have ii = and hence G[A, w = 0] = ef ^W^W = 1 
for any A. Thus (J6j> is correctly normalized. 

For the more rigorously minded reader, another deriva- 
tion of (|4]) and ([5]) is presented in appendix [A] It avoids 
the use of path integrals with unclear convergence prop- 
erties and takes into account the initial condition. 



III. EXAMPLES 
A. Stationary velocity distribution and propagator 

As a first application, let us re-derive the well-known 
probability distribution for the velocity in the case of 
stationary driving, w(t) = vt. 

To obtain the generating function of the velocity distri- 
bution at tp, we set X(t) = XS(t — t ) in ©. The solution 
of © is HI 

As already derived in [22|, for w(t) = v one gets 



l(t)w(t) = -«ln(l - A), 



(15) 



and hence G(A) = (1 — A) This generating function 
yields the probability distribution 



P{u) 



T(v) 



(16) 



which is the well-known result for the stationary velocity 
distribution 0, Q . 

Using the same method, we can obtain the 2-time ve- 
locity probability distribution. For X(t) = X\S(t — t\) + 
X2&{t — t^), with t\ < t2, the solution of ([5]) is 



u(t) 



1+- 



1— "J. "J- !£ — "i"- —a- - e t!- 

(1 + M)*2= 1 +Ai(l-A 2 )e*2 



t > t 2 
h <t<t 2 

t < tl. 



(17) 



As already derived in [22|], for w(t) = v one gets 
u(t)w(t) = 

~! _ e -(t a -ti) 



-v In |l — Ai — A2 + Ax A2 
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1 — Ai — A2 + A1A2 



1 



,-(t 2 -ti 



and using (|6]) 

G(A 1 ,A 2 ) = { 

JJ (18) 

Taking the inverse Laplace transform, we obtain the 2- 
time velocity distribution 



P(ui,U 2 ) 



2I»sinh§" 



-l+v 



\Ju\u-i 
sinh ^ 



where u\ :— u(t\), ii2 '■= iife), T t-x — t\ > and I a 
is the modified Bcssel function. This formula generalizes 
the quasi-static result of 22] to arbitrary v. Dividing by 
the 1-point distribution P{u\) given in (|16[) , one obtains 
a closed formula for the ABBM propagator for velocity 
v > 0: 



— l+v 



P(«a|«i) 



2 sinh 



—I-i 



v \ sinh Z J 



+ «2 



(19) 

Using this result and the Markov property of equation 
(O, n-point correlation functions of the velocity can be 
expressed in closed form as products of Bcssel functions. 



B. Velocity distribution after a quench in the 
driving speed 

Now let us consider a non-stationary situation. As- 
sume that the domain wall is driven with a constant ve- 
locity v\ > for t < 0, which is changed to V2 > for 
t > 0. One expects that the velocity distribution interpo- 
lates between the stationary distribution for vi at t = 
and the stationary distribution for V2 for t — > 00. In this 
section, wc will compute its exact form for all times. 

For the one-time velocity distribution, A(t) = A5(t— to) 
and the solution of §5§ is unchanged, given by (|T4")) . 

Now, using w(t) = {v\ — V2)0{—t) + v 2 , one gets 



.(t)w(t) 



viXdt 



A + (1 - A)e-(*"*o) 



V2X dt 



(vi - v 2 ) In 1 + 



1-A 



A + (1 - A)e-C*-*o) 
v 2 ln(l-A). 



Thus, with the help of 



G(A) - eMM = [1 _ A (l - e-*")]" 1 "" 2 (1 - A)" 
Inverting the Laplace transform, one obtains 

J (l - e -*0)«l-«2 



(20) 



P(u(t )) 



X 1F1 I V 2 - Vl,«2, 



1 - e*° 



(21) 




FIG. 1. Decay of velocity distribution after driving was 
stopped. Curves (from right to left): Results for to = 
0,0.4,1,2 from (|22[) . Bar charts: Corresponding simulation 
results, averaged over 10 4 trajectories. Initial driving velocity 
(for t < 0) was v = 10. 



An interesting special case is when the driving is turned 
off at t = 0, i.e. v\ = v and v 2 = 0. According to 
(jlip . the particle will continue to move forward until it 
encounters the first zero of u — F(u) — m 2 [u — w(0)]. 
Correspondingly, we expect that the velocity distribution 
decays from the stationary probability distribution at t < 
to a (5-distribution at zero at t — > 00. The explicit 
calculation for ii :— u(to) yields: 



P(u) = (l-e- t0 ) v S{u) 



(22) 



1 - e to 



The S(u) term gives the probability that the motion has 
stopped at time to, 



P(u(t ) = 0) = (1 - e-^y 



(23) 



As expected, this is zero at to = and tends to 1 as to — > 
00. Correspondingly, the distribution for the relaxation 
time T, i.e. the time for the particle to stop moving from 
the stationary driving state at velocity v, is given by 



P{T) - k 



P{u(t a ) = 0) = ve- T {l - e-' 1 ') 



-T\-l+v 



t =T 



The term in (|22p not proportional to the (5-function (once 
normalized) gives the conditional distribution of veloci- 
ties assuming the particle is still moving. Its form com- 
pares well to simulations, see figure [TJ 

Using (|20[) . one also sees that the mean velocity inter- 
polates exponentially between the old and the new value 
of the driving speed, 



tk = <9a| a=0 G(A) = i> 2 + (ui - v 2 )e~ 



(24) 



These results are valuable since they provide a tool to test 
the validity of the ABBM model in different experimental 
protocols. In application to Barkhausen noise, one could 
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perform experiments where the driving by the external 
magnetic field is stopped at some time. This would allow 
to verify e.g. (|22j) experimentally, since the velocity in our 
model is the induced voltage in a Barkhausen experiment. 
This would be one of the first checks on whether the good 
agreement between the ABBM theory and experiments 
persists in the non-stationary case. 



C. Non-stationary avalanches 

Using similar techniques, one can treat the case of a 
finite jump from to w in the location of the confin- 
ing harmonic well in ([T]), w(t) = w6(t) equivalent to a 
"kick" w(t) — wS(t). For t < the particle is at rest, 
and the quench at t — triggers exactly one avalanche. 
Its size is given by S = L u(t)dt and its duration T by 
the first time when u(T) = 0. Note that this avalanche 
occurs as the non-stationary response to a kick of arbi- 
trary size, a problem a priori different from the station- 
ary avalanches studied previously 0, H2, HH for small 
constant drive w(t) = v = + . In this section, we will 
derive the distribution of avalanche sizes and durations 
for arbitrary kick sizes w. 



1. Preparation of the initial condition 

The assumption G[X,w = 0] = 1 which we made in 
section Hi Bl implies that the initial condition at ti, which 
is the lower limit of all time integrals in the action and in 
(JBJ), is u{ti) = 0. This means that the particle is exactly 
at rest for t > ti if w(t) = for t > ti. Furthermore 
to assure that the particle will not revisit part of the 
trajectory, we demand u(t) < u(ti) for all t < ti. One 
protocol with which this can be enforced is: Start at 
some time t\ <C U at an arbitrary position u{t\) <C 0, 
and take w(t) = for all t g [tx,U]. Then u(0) will 
be almost surely positive. Thus, between t\ and i,, the 
particle will move forward until it reaches the smallest u 
where F(u) — m 2 u = 0. Since t\ -C ti, almost surely it 
will reach this point before ti and thus be at rest at ti. 
This choice of initial condition is equivalent to choosing 
a random configuration from the steady state for quasi- 
static driving at v = + . 



2. Duration distribution 



Laplace inversion gives, denoting u := u(to), 



P{u) 



e'O -1 



S(u) 



1 



2 sinh a 



sinh ^ 



The mean velocity 



u(t ) = d\\ x=Q G(\) = we' 



(26) 



(27) 



decays in the same way as in (|24[) for stopped driving. 
However the probability distributions of it (to) arc differ- 
ent, as can be seen by comparing (|26|) and (f22|) . The 
probability that ii(to) = 0, i.e. that the avalanche has 
terminated at time T < to, is obtained by taking the 
limit A — > -co in ([25]), which gives the 8- function piece 
in (1261). 



(28) 



P(ii(t ) = 0) = P(T < t a ) = exp I 



Note that this procedure requires P(u < 0) = which is 
the case here. 

Correspondingly, the probability density for the 
avalanche duration T is given by 



PiT) ^ S5 



P(i(io) = 0) 



w exp 



(2 sinh f) 



(29) 



We observe that for infinitesimally small quenches w, one 
recovers - up to a normalization factor - the distribution 
obtained in [22I [23| for avalanches at stationary, quasi- 
static driving, with the universal power law T~ 2 for small 
times Q: 



p(T) := d w \ w=Q P(T) 



1 



(2 sinh f; 



(30) 



Hence, the non-stationary character is not important in 
that limit. 

For finite w > 0, the mean avalanche duration is ob- 
tained from 



log w. 



T(w) =1e- e w ~Ei(-w) + log(to) 



It behaves as T(w) ~ w\n(l/w) at small w and diverges 
logarithmically for large w. In the latter limit, the distri- 
bution of T := T—hiw approaches a Gumbel distribution 



P(f) 



-T - 

e e 



on the interval Tg [—00, 00] , as if the duration were given 
by the maximum of w independent random variables. 



First, let us derive the exact distribution of avalanche 
durations following a kick. The generating function for 
P(u(to)) at time to > is obtained as in the previous 
section as 

""^""( uc'-y) < 25) 



3. Joint size and duration distribution 

One can now proceed to a more general case, and com- 
pute the joint distribution of avalanche durations and 
sizes. We again calculate the generating function 

G(Ai,A 2 ) = e A i s + A ^(*°). 
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where S :— J °° u(t)dt is the avalanche size. The solution 
of (|5|) for X(t) = Ai + A 2 <5(t — t ) is given by 

m = \{±- v / i-4A 1 ) + 

e^ 1 - 4Al < t - to Vl-4AiA 2 0(*o - 1) 
+ VI - 4Ai - A 2 [l - e v / T rr 4Ar(*-*o)] ' 

Since the driving is w(t) = wS(t), we obtain from §6§: 
G(\ 1 ,\ 2 )=e wZ ( Xl ' X2 \ (31) 



Z(A 1 ,A 2 ) = -(1-Vl-4A 1 ) 



e -^ 1 - 4Al *Vl-4AiA 2 



(32) 



VI -4Ai - A 2 (l - e -v / T rr 4A7to) ' 

For A 2 = this gives the distribution of avalanche sizes 
S for arbitrary kick size w, 

w 2 S w s 



P{S) = 



■ exp 



45 



(33) 



As it should, this coincides with the distribution obtained 
for quasi-static driving, v = + [44| . 

In the case of a non-stationary kick, we can obtain 
more information on the avalanche dynamics by con- 
sidering the joint distribution of avalanche sizes S and 
durations T. As above, the probability that u(to) = 
and hence the probability that the duration T of the 
avalanche lies in the interval ]0;to] is given by the limit 
A 2 — > — oo. Thus, the joint probability density P(S,T) 
of sizes S and durations T satisfies 



to 



dS / dTe Mb P(S,T) 
Jo 

exp I -(1 - v 7 ! -4Ai) - W y 



2 v v 1 - e -Vl-4Ai* 

Deriving with respect to to, we obtain 

/>oo 

/ dSe xs P(S,T) = 
Jo 

w(l - 4A)e^( 1 ~ v/T3IIcoth i VT^II) 
(2sinh|Vl —r 4A) 2 

which for A = reproduces (|2"9"j) . This implies the scaling 
form EH: 



(34) 



P(S,T) = e--f(S/T 2 ) 



lt: 



2iL — ^-\/scoth\/s 

we 2 se t 1 v v 



(35) 
(36) 



T 4 (sinhV^) 

Although no formula to invert the Laplace transform in 
a closed form is evident, one can, for example, calculate 
the mean avalanche size for a fixed value of the avalanche 
duration, 



S(T) = 



f™dSSP{S,T) 
J™dSP(S,T) 

4 - wT - 4 cosh T+(2T + w) sinh T 
coshT- 1 ' 



(37) 



As w 0, this has a well-defined limit 

S(T) = 2Tcoth^ -4. (38) 

Eg . (|38p reproduces the expected scaling behaviour 0, 
|21| . S(T) ~ T 2 for small avalanches. This is apparent 
in (|35|) . since the e _: r factor can be neglected for small 
S. The new result in Eq. (|3"8"|) predicts the deviations of 
large avalanches from this scaling, and shows that they 
obey S ~ T instead. This is in qualitative agreement 
with experimental observations on Barkhausen noise in 
polycristalline FeSi materials 0, [n], [2^|. It would be 
interesting to test quantitative agreement of (|38|) with 
experiments, as well. 

We can also obtain the large-T behaviour at fixed S 
(fixed A) since in that limit 



/>oo 

/ dSe xs P(S,T) w w(l 
Jo 



4A)e' 



(39) 



This implies 



u>(2T + w) [(2T + w) 2 - 65] 4« -f 

2V^5 7 /2 



P(S,T) 

(40) 

Note that (|39p is also valid at fixed T and large negative 
A, hence (|40p also gives the behaviour for S -C T 2 at fixed 
T. One notes some resemblance with (|33|) . 

We now consider the limit of a small kick w — > 0. 
Eq. ((HU gives 

P(5,T) = W p(5,T) + 0( w 2 ), (41) 

where p(5, T) can be interpreted as an avalanche size and 
duration "density", satisfying 

(1 - 4A) 



POO 

/ dSe xs p(S,T) 
Jo 



(2sinhfVl-4A)' 
This Laplace transform can be inverted: 



p{S,T) = e 
g(x) = LT^ 

+ 00 



- S/4 ^.9(5/T 2 ) 



(sinh Vs) 



(42) 

(43) 
(44) 



h(x) = ^ (1 - 27r 2 n 2 a;)e~ 



E 



2m 2 e 



,3/2 



n— — oo 



We have used 



l/(sinhVs) 2 - Note 



-oo (s+n 2 7r 2 ) 2 

that p{S,T), as a size density, is normalized to 
J °°dSp(S,T) = p(T), given in flU, since a fixed du- 
ration T acts as small avalanche-size cutoff. The total 
size density p(S) = J dT p(S,T) = 2 ^ s3/2 exp(-|) is 
not normalized, since w, which acts as a small-scale cut- 
off in (|3"3")l , has been set to 0. 

Finally, note that (f34|) allows one to go further and 
compute any moment as well as, by numerical Laplace in- 
version, the full joint distribution P(S, T). This is shown 
in figure HJ 
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FIG. 2. Joint density p(S, T) of avalanche sizes S and dura- 
tions T in the ABBM model, obtained by numerical Laplace 
inversion of (|42I43I44|) . The red line is the mean size S(T) for 
a fixed duration T given in (|38[) . 



4- Avalanche shape following a pulse 

We consider now the joint probability of velocities at 
two times < t\ < t 2 following a pulse at time t = 0. By 
©, its generating function is 



3 Ai«(ti)+A2'!i(t2 



Llfi(0) 



where u(0) is the 2-time solution (|17[) . We are interested 
in P(u(ii), w^) = 0) obtained by taking A2 — > —00: 



duie AlUl P('iti,0) = exp 



Aie t i+(l-A 1 )e«2 u 
(l + AiJe'i-Aie'a c 



1 



We use that LT7i„e d+ ^ 



with d 
1 



= a = ^/(e 41 - 1)2 ail d b = 1J ^ T - 1 + 

Taking dt 2 and setting ti = T we find the joint 



l-e-'i 

probability distribution of the avalanche duration T and 
the velocity iiiti) = u\, 



P{u\,T) = -<9 f2 fee d v /^/i(2y^")e" 
1 V wiii ^ 

T-ti^l 2 2sinh S 1 



[2sinh(^-)] 



\t 2 =T 

\Jwiii 
sinh ^ 



(45) 

Dividing by P(T) given in (|29|) . we find the conditional 
probability for the velocity distribution at t± for an 
avalanche of duration T . In particular, we get the av- 
erage avalanche shape, 

4sinh(^-)sinh(2^*i) 



u{ti) T 



sinh^) 



sinh(I^i) 
sinh(^) 



(46) 



For w — > one recovers the stationary avalanche shape 
obtained in [22|, [24| . On the other hand, avalanches fol- 
lowing a pulse of size w > have an asymmetric shape, 
since u(t = + ) = w. This should provide an elegant 
way to discriminate between the two situations experi- 
mentally. 



D. Power spectral density and distribution of 
Fourier modes 



In signal analysis, an important observable used to 
characterize a time series is the power spectral density 
P(uj) defined as 



P(w) 



lira — 

T->oo T 



,772 






/ pilot 


u(t) - u(t) 


At 


J-T/2 







(47) 



This gives a measure for the abundance of the frequency 
component u> in the time series u(t). For a stationary 
signal where the 2-time velocity correlation function only 
depends on the time difference, (|47[) is equal to its Fourier 
transform: 



e mt u(0)ii(i) dt. 



(48) 



For driving with constant velocity w(t) = vt, one knows 
l46j ii(0)u(t) = ue - '*' and hence the power spectrum 
for the velocity in the ABBM model is 

2v 



1 



(49) 



We can now proceed further and obtain the probability 
density of each Fourier component. We consider ([6]) with 
X(t) = X cos tut 6 (T — t) 9 (t) where T is a large-time cutoff. 
To solve (O with this choice of A, we substitute u(t) = 



giving Mathieu's equation, 



4>"{t) - Q - AcoswtJ <p{t) = 0. 

This is to be solved with the boundary condition u(T) = 
0, i.e. <t>'{T) = -\<t>{T). _ 

The general solution is a linear combination of two 
Floquet solutions 



m=^Pi{t) + e-^P 2 {t), 



(50) 



where P\ t 2{t) are periodic functions, pi = /i(A,w) is re- 
lated to the conventionally defined Mathieu characteristic 
exponent v{a,q) (in the notation of 47j) by 



1 



2A 
7a 



When A is real and close to 0, pi is real, has the same sign 
as A, and is odd in A. Thus, for < t <C T, the solu- 
tion (j)(t) given in (|50|) is dominated by the exponentially 
decaying term 

4>{t) « e-^^Pit), 
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with Pit) = P\^{t), depending on the sign of A. Thus, 
for < t <C T we have 



1 P'(t) 
u(t) = --rt\M,w) + 



P(t) 



(51) 



In order to evaluate one needs to integrate ii(t) over 
t from to T. Since P(t) is periodic, its contribution 
vanishes for each period 



u(t) dt = — 

10 



i-A*(|A|,^ 



< s < T. 



(52) 

For constant driving, w(t) = ui and T ^> ^j, one thus 
obtains using © 



3 A/ T ii(t)coswt _ vf^u(t) 



L )]+o(T) 



. (53) 

As expected by symmetry, this is an even function in 
A. It remains real as long as the Mathieu exponent v 
is purely imaginary, which is the case for |A| < X c (lo). 
One can interpret the corresponding Mathieu functions 
as Schrodinger wavefunctions in the periodic potential 



V(x) 



1 



A cos(wa;) 



The region |A| < X c (lu) is the region where the energy 
E = is outside the energy band(s) of this potential, 
and all wave-functions are evanescent. At A = ±A C , 
one has v = and for |A| > A c , i.e. outside the "band 
gap" , the expectation value on the left-hand side of ([53]) 
does not exist. This indicates that the distribution of 
J u(t)cosiot has exponential tails for any to > 0. The 
exponent of this tail can be computed in terms of the 
so-called Mathieu characteristic values a r and b r |47j . 

Furthermore, from (|53"j) one observes the scaling be- 
haviour of the cumulants 



u(t) cos cot 



T, 



(54) 



which reminds of the central limit theorem. 

Taking two derivatives of (|53[) with respect to A, and 
using dq\ _ Q+ v(— b 2 , q) = 2 b(b 2 \i) > one ver ifi es once more 
(|4U|) . However, (|53|) goes beyond that and gives the full 
probability distribution of each frequency component of 
the time series u(t). 

With this, we conclude our examples on the "classical" 
ABBM model and move to generalizations which can be 
treated by our method, as well. 



IV. ABBM MODEL WITH SPATIAL DEGREES 
OF FREEDOM 

An interesting generalization of the ABBM model (p} 
is a model with spatial degrees of freedom, (e.g. an ex- 
tended clastic interface in dimension d > 0), but subject 



to the same kind of disorder as in the ABBM model, i.e. a 
pinning force correlated as a random walk. 

An interface was studied in [22J for quasi-static driving 
and it was found that the global motion (i.e. the motion 
of the center-of-mass of the interface) is unchanged by the 
elastic interaction. An instanton equation for the other 
Fourier modes was derived, but solving it remained a 
challenge. 

Here we extend these results to arbitrary driving veloc- 
ity. We first study the simpler case of only two elastically 
coupled particles, and present a direct argument to show 
that the center of mass is not affected by the elastic inter- 
action and is the same as for a single particle, i.e. model 
([T|) in a rescaled disorder. For two particles the instan- 
ton equation is simpler and more amenable to analytic 
studies, which allows us to see how local properties (such 
as the velocity distribution of a single particle) are mod- 
ified. In the last part we come back to the interface and 
show a non-renormalization property of the theory valid 
for any driving velocity. 



A. Two elastically coupled particles in an 
ABBM-like pinning-force field 

The model we analyze in this section is a 2-particle 
version of ((T|): 

d tUl (t) = Fi(«i(i)) - m 2 [ui{t) - w(t)] 

+k[u 2 {t)-u 1 (t)}, 
d t u 2 (t) = F 2 (u 2 (t)) - m 2 [u 2 (t) - w(t)] 

+k[ Ul (t)-u 2 {t)}. (55) 

We assume Fi(ui), F 2 {u 2 ) to be independent Gaussian 
processes with correlations as in ©, i.e. 



[Fi(u) - Fi(u')] = [F 2 (u) - F 2 {u')Y = 2a\u - u'\. 



1. Center-of-mass motion 

From (f5"5"j) . we obtain the equation of motion for the 
center-of-mass velocity s(t) = h[ui(t) + ii2(t)]: 

d t s(t) = \d t \F x {u x {t)) + F 2 (u 2 (t))} - m 2 [s(t) - w(t)} . 

(56) 

To better understand the effective noise term 
dt[Fi(ui(t)) + F 2 (u 2 (t))}, let us compute its gener- 
ating functional, 



G[A] = e f t Ht)idtlF 1 (u 1 )+F 2 (u 2 )] 

= e -ft,t> A(t)A(t')f9 t a t ,(|« 1 (t)-«i(t')|+|« 2 (i)-« 2 (i')l)_ 

Using monotonicity [3 |4(| of the trajectories we 
obtain 

G[A] = e^ A(t)2 T M*)+«2(i)] = e / t A(t) 2 fi(t)^ 
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Note that this is the same generating function as for a 
random pinning force F(s(t)) with correlations 



[F(s) - F(s')}= a\s - s' 



Thus, we can re- write (|56| as 



d t s(t) = d t F(s) - m 2 [s(t) - w(t)} 



(57) 



(58) 



with a rescaled disorder amplitude a' = § , reducing it to 
the same form as ((7]). 

This argument extends straightforwardly to any num- 
ber of elastically coupled particles, and to the continuum 
limit. Thus, we observe that the dynamics of the center 
of mass of an extended interface in a pinning-force field, 
which is correlated as a random walk, is equivalent to the 
1-particle ABBM model (P). 



2. Single-particle velocity distribution 

On the other hand, observables that can not be de- 
scribed solely in terms of the center of mass are more 
complicated. In order to obtain the joint distribution 
of the particle velocities Ui(t), ii2{t) one may follow the 
same route as in section Hi Bl We start from 

G[\x,\ 2 , W } = e / f Ai(*)ui(t)+A 2 (t)« 2 (t) 

= e ™- 2 I t [iil(t)+U2(t)]w(t) ^ j-gg-j 

where iii, u 2 are solutions of the coupled nonlinear dif- 
ferential equations 

- dtMt) + m 2 ui(t) + k [ux{t) - u 2 {t)\ 

-au^t) 2 = Ai(t), 

~d t u 2 {t) + m 2 u 2 (t) + k [u 2 (t) - ux(t)} 

-au 2 (t) 2 = X 2 (t). 

In contrast to (0, these can not be solved in a closed 
form even for simple choices of \\_ 2 . However, one can 
obtain a perturbative solution for small k around k = 0. 
To give a simple example, one obtains for monotonous 
driving w(t) = vt and one-time velocity measurements 
Ai,a(t) = Ai, a <5(t): 



G(Ai,A 2 ) = [(1 

l + ufc(Ai - A 2 ) 

+ 0{k 2 



Ai)(l-A 2 )p 
ln(l-Ai) 



ln(l - A 2 ) 



Ai(l-A 2 ) A 2 (l-Ai) 



(60) 



where we use rescaled units where k denotes k/m 2 in 
the original units. As one expects from the previous sec- 
tion, the correction of order k vanishes if one considers 
the center-of-mass motion, Ai = A 2 . If, on the other 
hand, one considers the 1-particle velocity distribution, 
i.e. takes A 2 = 0, one gets 



G(Ai,0) = (l-Ai)-" (1+fc) 



1 - vk- 



Ai 



1-Aj 



0(k 2 



(61) 
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FIG. 3. Single-particle velocity distribution P(ii) in the 2- 
particle toy model for weak elasticity. Histogram: Numerical 
results from simulations for k — 0.2. Dashed line: Stationary 
distribution in absence of elastic coupling {k = 0). Solid (red) 
line: O(k) result from (p2l) . 



The Laplace transform can be inverted, giving 

e-^ Ul - 1+v 



P{ui) 



T(v) 



x 



(62) 



x jl + k [v — iii + v In iii - vi'{v)} +0(k 2 )} 



where ip(x) = j^j is the digamma function. Simulations 

for small k confirm this result (see figured]). The next or- 
der in k can likewise be calculated, however the resulting 
expressions are complicated and not very enlightening. 

A non-trivial consequence of (|62[) is that the power-law 
exponent of the distribution P(iii) for small velocities 
changes from ii~ 1+v to { l - 1 + v ( 1 + k ) . 



B. Continuum limit and non-renormalization 
property 

Let us now consider a <i-dimensional interface in a d+1- 
dimensional medium with a generic elastic kernel g xy> 
such that in Fourier = m 2 . Local elasticity corre- 
sponds to g^ 1 = q 2 + m 2 . The corresponding generaliza- 



tion of ([T]) is 



d t u xt = F(u xt ,x) 



9xy u Vt + ^xt- 



(63) 



For the remainder of this section, we write function 
arguments as subscripts in order to simplify notations 
(i.e. u xt := u(x,t)). The source X x t > for the field ii 
is a positive driving, and is related to the velocity of the 
center of the quadratic well w by X x t = g xx ,w x 't- 

The pinning force is chosen gaussian and uncorrelatcd 



m x, 



F(u, x)F(u>, x') = 5 d (x - x')A(u, u'). 



(64) 



In the u direction, analogously to we assume 

Brownian correlations, i.e. uncorrelated increments: 
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d u d u 'A(u, u') = S(u — u'). This does not fix F uniquely, 
with e.g. two possible explicit choices in ([71)) and (|75|) 
below. However, differences only arise for the position u 
but not for the velocity u, as will be discussed below. 

Let us write the MSR partition sum in presence of 
sources, 



We can now write down the effective action T[u, u] 



G[X, A] = / V[u\V[u]e- s[ ^ + ^ Xxt<llct+ ^ 



The generalization of the MSR action (|T3|) to this situa- 
tion is 



S[i 



u xt I d t ii xt + I g X yU yt - (Tu xt u xt 

xt J y 



(65) 



To arrive at (|65[) we have again assumed forward-only 
trajectories u xt > 0, guaranteed if X x t > and u xti > 
at some large negative initial time ti. 

The solution in section III Bl generalizes straightfor- 
wardly to 



G[A, w] = eft = e-k W A ** , (66) 

where [A] is defined as the solution of 

dtuil - I g^4 3 t + = ~^t- (67) 

Jy 

In principle, this can be used to compute any observable 
of the d-dimcnsional theory. In practice, the equation 
(f67|) for u is hard to solve analytically for most cases. 

In the remainder of this section, instead of discussing 
specific examples, we show a conceptual consequence of 
(|66j) : The action (|65|) does not renormalize. The effec- 
tive action T is equal to the microscopic action S in any 
dimension d. 

According to (|66|) . the generating functional for con- 
nected graphs W[X, A] evaluates to 

W[X,X]=lnG[X,X] = f u£>[X]X xt . 

Jxt 

To perform the Legendre transform from W to the effec- 
tive action T we introduce new fields ii xt [X, A], and 
u xt [X, A], defined by 



u x t 



SW[X, A] _ w 

= = U Tt 

SX xt 
SW[X, A] 

sx xt 



Su { x %[X}~ 



X 't' 3X xt 



X x 't'- 



(68) 
(69) 



Here and below we drop the functional dependence on 
the sources when no ambiguity arises. Eq. (|68|) shows 

that u X f [A] is really the field u xt appearing in the effec- 
tive action, hence (|57j) allows to express the field X x t (on 
which W depends) in terms of u xt (on which T depends). 



T[u,u] = 

Jxt 



u x tX x t 



UxtX x t 



w 



u xt X xt since W = u xt X xt 

xt J xt 

Uxt (d t Uxt ~ J 9xyUyt + Ou\ t 



u xt I d t u xt + I g X yU yt - o-u xt u xt 

xt \ J y 

S[u, u]. 



(70) 



This is exactly the same as the bare action S in (|65p . This 
non-renormalization of the action for the particle veloc- 
ity in ABBM-like disorder is also consistent with a 1-loop 
calculation using functional RG methods (see appendix 
|B|) . It is a very non-trivial statement, and shows that, in 
some sense, the MSR field theory for monotonous mo- 
tion in ABBM-like disorder is exactly solvable in any di- 
mension. The monotonicity assumption implies that the 
derivatives arising in the formulae above must be per- 
formed in the neighborhood of a strictly positive driving 
source X x t > 0. Using the relationship 



. r , r, Uxt exp T I ', ., Xx't'Ux'V 

u x t[X,X] = == = i 

eX P L/x'f X *'t'Ux't>\ 



(71) 



(where the average is performed in presence of w = A) 
one sees that (|69p maps positive A onto positive ii. On 
the other hand, the condition A > can be expressed 
using X x - 



t)it.r 



as 



Uxt < 



d t u x 



+ Jy9xyU 



2(TU, 



(72) 



We conclude that the effective action T[u, u] is given by 
the bare action S in the sector of the theory where it > 
and (|72|) holds as a necessary condition. In no way this 
implies that T = S for values of the fields where this 
monotonicity assumption does not hold. The case of non- 
monotonous motion and/or non-monotonous driving is 
highly non-trivial and will be studied elsewhere. 

In the following section, we shall see how this result 
generalizes to the field theory of the position u(t), where 
the relationship between S and T is slightly more com- 
plicated. 



V. FIELD THEORY FOR THE POSITION 
VARIABLE 

So far, we have considered observables that can be ex- 
pressed in terms of the ABBM velocity ii(t), or in case of 
a manifold u{x,t). Here we consider the position u(x,t) 
itself. One can then formulate the MSR path integral in 
terms of u and it, analogous to (|12[) . This was done for a 
<i-dimensional interface in short-ranged disorder in [23| . 
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as a starting point for a d c — d-expansion. Here we focus 
on the simpler and solvable case of the ABBM model, 
where the MSR path integral reads 



S[u, u]= u(t) [d t u(t) + m 2 (u{t) - w(t))] 



A(u(t) 7 u(t'))u(t)u(t'). 



(73) 



Here, A(u,u') = F(u)F(u') is the disorder correlation 
function. One mathematically simple choice is to assume 
the random force F(u) to be a one-sided Brownian mo- 
tion and restrict to u > 0: 

A(u,u') = 2crmin(u,u') = a(u + v! - \u - u'\). (74) 

Another common choice is the two-sided version, i.e. a 
Brownian motion on the full real u axis pinned at F(u = 
0) = 0. With either choice, however, the random force 
is non-stationary and one loses statistical translation in- 
variancc. This is unnatural for certain applications, for 
example approximating extended elastic interfaces above 
the critical dimension. In this context, one chooses a 
stationary variant of (f7"4")) . 



A(u,u') = A(u-u') = A(0) -a\u-u'\. 



(75) 



Since a stochastic process F(u) can only satisfy (|75[) for 
all u in some limit, we always assume (|75p to be regular- 
ized at large \u — u'\. 

For observables that can be expressed in terms of the 
velocity ii, only dtdt>A(u(t), u(t')) enters the MSR action 
(cf. section Hi Bp . Hence, choosing (|74[) or (|75|) yields the 
same result (flU)) . However, the choice does matter if one 
is interested in observables depending on the position, 
like the mean pinning force f p := m 2 [u(t) — w(t)}. 

In contrast to the velocity theory discussed in previ- 
ous sections, fixing a distribution of positions u(tj) as 
the initial condition is problematic. Indeed, in general 
one cannot exclude that this initial condition leads to 
backward motion ii{ti) < for some realizations of the 
disorder. Hence for the stationary Brownian landscape 
(|75j) we will choose t, = — oo and assume that the driving 
w(t) > is such that at fixed times the initial condition is 
forgotten, as discussed in section HlI C II We claim that 
then 



(76) 



1 - "I / A W 



G[X,w] = e& A ('M*) 

= e m 2 f t Ht)w(t)+f™[f t \(t)} 2 



where all time integrals are over ] — oo, oo[. The function 
u(t) = —dtu(t) where u(t) is solution of 



d t u(t) — m u(t) + au(t) — au(t)u{— oo) 



t'>t 



A(f). 
(77) 



In the particular case of the one-sided Brownian land- 
scape (|74p we only consider the initial condition u(ti) = 
0. Since F(0) — in that case, for w(U) > and w(t) > 
the motion will be forward. Then the generating function 
G[A, w] in (|73j) takes a form analogous to (|6|) 



G[\,w] = e 1 ^ A(t)u(t) = e m S »'i " (tMt) , (78) 

where u(t) = —dtu(t) and u(t) is solution of (|S|). In 
the remainder of this section, we shall prove the above 
statements and then apply these formulae to determine 
the distribution of the single-time particle position u(t). 



A. Generating functional for stationary Brownian 
potential 

Using the assumption of monotonous motion, the dis- 
order term in the action (|73[) can be rewritten as 



A(u(t),u(t'))u(t)u(t') 



t.r 

A(0) 



u(t) 



u(t)u(t)u(t')sgn(t - t'). 



Following the same approach as in section lTlBl evaluating 
the path integral over u(t) in (|73p yields 



m 2 / f «(t) l0 (t) + ^l[/ t u(t)] 2 



V{u]e 



x 5 d t u{t) - m 2 u(t) + au{t) I u{i )sgn(f -t) + X(t) 



(79) 



Thus 



G[A, w] - A/e™ 2 /. Ht)Mt)+mf t , (80) 
where ii(t) is solution to the equation 

d t u(t)~m 2 u(t) + au(t) / u(t')sga(t' - 1) = -A(i). (81) 



Substituting u(t) := J°° u(t)dt, one recovers ((771) . 
u(— oo) is obtained from 



m / u(t) = —m u(— oo) 



X(t')dt'. (82) 



Note that u(— oo) vanishes for A such that J., X(f) = 0. 
These are exactly those observables which can be ex- 
pressed in terms of the velocity (or, equivalently, position 
differences) . 

As in section mi JV in (|80p is the normalization of the 
path integral and the Jacobian of the operator inside the 
(^-functional in (|79p . It is independent of w(t), but we 
cannot fix its value at w(t) — const as we did for the 
velocity theory in section [TT| Even if one keeps w = const 
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for a long time, the distribution of u will remain non- 
trivial (unlike the distribution of ii, which will become 
8(u)). Here, to fix Af we compare to the disorder-free 
solution (ct = 0) for which the trajectory u{t) is deter- 
ministic and satisfies (|80[) with Af = 1. Hence, we can 
write Af as a ratio of functional determinants arising from 
the ^-functional, 



cedurc as in section llVBl one obtains the effective action 



T[u, u] 



hi 



A(0) 



d t u(t) + m 2 u(t) - au(t) j u(t')sgn(t' - t) 

2 



(t) 



- In 



i - — 2 ! m 

m 1 J t 



Af~ 1 = 



det(d t -m 2 - S T ) 



det(9 t — m 2 ) 
Here, R is the disorder-free propagator 



= det(l + RE). (83) = S[u, u] - In 



R := {d t 



Rtx,t 2 — 8(ti — t 2 )e 



(84) 

and E is the disorder "interaction" term, or " self-energy" 



<?S(h - t 2 ) j ' u(t')sgn(t a - t') 
+au(t 2 )sgn(t 2 -t 1 ). (85) 



By explicit computation (see appendix [Bj , one verifies 
that 



tr (RE) r 



~ ( u{t) 
m 2 J t 



and hence 



det(l + RT,) = exptrln(l + RT,) = 1 - 



u(t) 



From (|81|) . one further knows that J t u(t) = -^j- L \(t). 

In total, this proves the expression (|T6[) for the station- 
ary case, 



G[X,w] 



1 



m 4 .,t 



X(t) 



(86) 

One sees again that for observables expressed in terms 
of the velocity, where J t X(t) = 0, the simpler expression 
(HU) is recovered. 

In the language of pcrturbativc field theory, the non- 
trivial functional determinant signifies non-vanishing 1- 
loop diagrams [5l|. This is in contrast to the theory 
for the velocity (section [II B[) . where all observables were 
given by tree-level diagrams. These loop corrections 
mean that the non-rcnormalization property discussed 
in section IIVBI has to be amended when considering the 
particle position in a stationary potential. After renam- 
ing the driving w to A = m 2 w, the source for the field u, 
the generating functional for connected correlation func- 
tions becomes 



W[X,X] = u t [X]X t + 



A(0) 
2m 4 



-In 1 



X, 



where Ut[X] is solution of ([8"Tj) . Following the same pro- 



1 - —2 [ «(*) 



(87) 



We thus see that the property T = S seen for the velocity 
theory is only changed by a simple contribution from the 

1- loop corrections. The equal-time part of the u n term 
of these loop corrections coincides with a previous result 
in [52|. 

In fact, this calculation can be extended to the d- 
dimensional interface with elastic kernel g q of section 
IIVBI There too it ensures that for the position theory, 
and monotonous driving, T differs from S only via the 
logarithm of a (one-loop) functional determinant. Thus, 

2- and higher-loop corrections to correlation functions 
and the effective action vanish. Its expression is partic- 
ularly simple in the case of a uniform X x t = X(t) leading 
to a uniform saddle point u xt = u(t): 



loop 



d d q 
(2^p 



In 



m Jt 



T-S\ 



uniform 



- — T a 



d d q 
{2n) d 



In 



1 - ag q j u(t) 



L d is the volume of the system. Details and a more gen- 
eral discussion are given in appendix [B] appendix [C] and 



B. One-sided Brownian potential 

It is instructive to give for comparison the solution 
for the simpler case of the correlator (|74|) . Using the 
assumption of monotonous motion, the disorder term in 
the action (|75|) can be rewritten as 

a I n\in (u{t),u(t'))u{t)u{t') = 2a u(t)u(t) / u(t') 
Jt.t' Jt Jt'>t 

(89) 

Following the same approach as in section lHBI evaluating 
the path integral over u(t) with initial condition u(ti) = 
in (f73|) yields equation (|78|) . where u(t) is solution to the 
equation 

d t u{t) - m 2 u(t) + 2au{t) I u{t') = -X(t). (90) 

Jt'>t 



Note that as in section Hi B I the initial condition u(ti) = 
ensures that G[X,w = 0] = 1. Hence the functional 
determinant analogous to (|83[) is equal to 1 in this case. 
This is also checked by a direct calculation in Appendix 
151 For X(t) non-vanishing only around t ^> ti and w(t) 3> 
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w(ti), we expect that the influence of the initial condition 
is negligible. In this particular limit, ([75} should hold 
independently of the initial condition. 

Introducing u(t) := f t , >t u(t'), ([501) gives the following 
equation for u(t): 

d t u{t) - m 2 u{t) + au(t) 2 = - [ X(t'), (91) 

Jt'>t 

where we used that u(t) — > for t — > +oo (we recall that 
ii(t) must vanish at both ±oo). 



C. Example: Single-time position distribution 

To give a simple application of ([76]), we compute the 
distribution of the position u(t) at a single time. To do 
this, set X(t) = X8(t — t ) in ([731 . For the Brownian case, 
one obtains 



u(t) 



A(l - AX)9(t a - t) 



sinh 



VT^4A(t-to) 
2 



2 



- 4 A cosh 
For the stationary case (|77p . u(t) reads 

A (1 - A) 2 e-^-*»)( 1 - A >6»(t - *) 



u(t) 



[ e -(*-*o)(i-A) - A]' 



In both cases, the 6* functions come from causality, since 
the driving w(t) for t > to cannot influence the measured 
position u(to). Hence both u(t) and u(t) = —dtu(t) must 
both be identically zero for t > t . 

Let us assume a constant driving velocity, and write 
w(t) = v(t — ti) + Wj. Then, for the one-sided Brownian 
with u(ti) = and u>i > we have 

G(A) = e^o" = e^i di ' fi (*')«*'_ 

This leads to a complicated formula which simplifies in 
the limit i, — > — oo at fixed w{to), 



G(X) = 



-2A 



^oi(l_VT=4A) 



, 1 - 4A - V1-4A 
For the stationary case (restoring units), this is 

'■J t , u{t')vt' dt' + |WA 2 



(92) 



G(X) 



m 4 



1 



-A 



A(0) ,2 



Inverting gives a valid distribution only for \ut — vto\ <C 
cr/A(0) which coincides with the cut-off which should 
be used to regularize the stationary Brownian landscape 

(USD- 



VI. GENERALIZATIONS 

In light of the interesting results obtained for ([T]), it 
is natural to ask whether our approach can be extended. 
In particular, one might want to replace the response 
function in ([1]) by a more general response kernel. For 
example, in order to model eddy currents which change 
the avalanche shape in real magnets [jl,[53|, one may want 
to include second-order derivatives in time. 

For this, it is useful to view the calculation in section 
III B I from another perspective. The equation ([5]) for u is 
identical to the saddle-point equation obtained from the 
action (| L3[) in presence of the source A by taking a func- 
tional derivative with respect to u(t). The result ([6]) is 
then the value of Z at the saddle point obtained by solv- 
ing ([5]) for the given choice of A. The other "coordinate" 
of the saddle point (which happens not to influence the 
value of Z in this case, however) is the field u(t), fixed by 
the equation obtained by a functional derivative of (|13l) 
with respect to u(t), 



d t u(t) + m 2 [u{t) - w(t)} - 2au(t)u(t) = 0. 



(93) 



This is the trajectory giving the dominant contribution 
to Z for a given choice of A. E.g. for X(t) = XS(t — to), 
u(t) is given by (|14[) ; for w(t) = vt the solution of 
converging to v at infinity then reads 



u(t) 



A 



1 - A 



-I*— *ol 



Note that it can also be obtained from the 2-timc 
generating function ([T5[) . e.g. for t > t as u(t) = 
c*A 2 lnG(Ai = A, A 2 )|A 2 =o,« 2 =«,ti=to' Indeed, since S = T 
for monotonous motion, the solution of (|93p identifies 
with (|71[) . i.e. the saddle-point approximation is exact. 
We thus see, as expected, that if we concentrate on small 
velocities (A — > — oo), the velocity on the dominant tra- 
jectory u(t) gets closer and closer to at to, but never 
becomes negative. 

Now, the action S generalizing ([TB"]) with an arbitrary 
response kernel R U ' is 



S[i 



(94) 



u(t) 



Rt t }u{t') - m 2 w{t) 



aii(t)u(tf 



The saddle-point equations read 



K 



-1,T~ 



u{t') - ou{ty - x{t) = o, 



Rt t }u(t') ~ 2<ru(t)u(t) - m 2 w(t) = 0. (95) 



For a general (bare) response function R, the last term 
in the action (|94[) is not exact, since we cannot assume 
monotonicity of each individual trajectory. However, as 
long as the saddle-point trajectory defined by (|95[) for 
some choice of A is monotonous (i.e. satisfies u(t) > for 
all t), it gives a well-defined approximation to the value 
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of Z for this particular A. Investigating the quality of 
this approximation is an interesting avenue for further 
research. 



VII. SUMMARY AND OUTLOOK 

In this paper, we have considered the ABBM model 
with a monotonous, but non-stationary driving force. Us- 
ing the Martin-Siggia-Rose formalism, we obtained the 
generating functional for the velocity from a field theory 
that can be solved exactly. This was illustrated on sev- 
eral paradigmatic examples (e.g. a quench in the driving 
velocity). Using our formalism, we also succinctly recov- 
ered previous results on the stationary case. 

An interesting direction for further research is trying 
to generalize these results to non-stationary dynamics 
of models which are not mean-field in nature, like d- 
dimensional elastic interfaces. Although some work has 
been done in that direction [54l - l57| . many questions re- 
main open. Another complication arises when adding 
non-linear terms to the equation of motion ([T]) or (|63j) . 
The effects of the KPZ term [Vu(x)] 2 have been dis- 
cussed in j58l - l60j . An anologous term but with a time- 
instead of a space derivative, i.e. a term u 2 , is related 
to dissipation of energy [6l[ and yields a toy model with 
velocity-dependent friction. This is important as a step 
towards realistic earthquake models, where it is known 
that instead of a constant friction coefficient one has a 
complicated rate-and-state friction law fl^4l4^ . For the 
hysteresis loop in the ABBM model, it would be interest- 
ing to extend our results to the case of non-monotonous 
driving. Unfortunately, this is not an easy task: We cru- 
cially used both the monotonicity of the particle velocity, 
u{t) > 0, and the one of the driving, w(t) > for sim- 
plifying the action and computing the path integral in 
section Hi Bl Without this assumption, neither the result 
^ nor the non-renormalization property in section TlVBI 
hold. Assuming the non-renormalization property, the 
mean velocity u(t) would be equal to its value in the 
system without disorder at all times. This can be seen, 
e.g. by taking d\ at A = in formula and using © 
and l|14p. However, in numerical simulations one observes 
that this property breaks down as soon as the driving is 
non-monotonous, hence at least the term proportional to 
u in the effective action is renormalized. We thus leave 
questions in this direction for future studies. 
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Appendix A: Derivation of the non-stationary 
solution in discretized time 

The path integral derivation of (j6]) in section |n] is, to 
some extent, formal and neglects subtleties like conver- 
gence issues and boundary conditions. To complement 
it, we provide here a rigorous first-principle derivation of 
© by discretizing the time axis. For a small time step 
St, we write ([7]) as follows: 

Uj+i - Uj = F (uj + Stv,j +1 ) — F(uj) 
St St 



+m (w j+ i - Uj+i) 
=> iij+i = X{uj+i) + km 2 Stwj+i + kiij, (Al) 
with fc _1 := 1 + m 2 St. 

X(v,j + i) := k{F (uj + Stiij + i) — F(itj)] is, by the 
Markov property of Brownian motion, a new Brown- 
ian motion with X(0) = and variance X(u)X(u') — 
2ak 2 Stmm(u,u'). Eq. (|Alj) is an implicit equation for 
iij+i, which has, in general, several solutions Uj+\ > 0. 
In fact, its solutions are the intersections of the Brownian 
motion X(iij+i) with the line km 2 8twj+i + kiij — iij+x- 
The true iij+i describing the motion of the particle is the 
smallest of these solutions. 

Hence, the conditional probability distribution for iij+i 
given uj is the first-passage distribution of Brownian mo- 
tion, given by 



P{u 3+ i\uj 



km 2 StWj-i-i + kiij 



o 2 



\/4:Trak 2 5tUj +1 



(A2) 

The Laplace transform of this expression, which is the 
conditional expectation value for e UUj+1 , is given by 

/>oo 



+ " 
g 2akSt 



^(l-Vl-4«CTfc 2 (5t) 



This can be rewritten as 



(A3) 



(A4) 



with u 



2akSt 

H) one obtains 



[l — VI — 4uak 2 5tj. Hence, iterating 



gE^Li \jU } St _ e m 2 J2f=i UjibjSt x gfiiito^ (A5) 



where itj is defined via the (backward) recursion 
u N+1 = (A6) 

_ i-yi-4(a j+1 ft + y^ <7<iv 

3 2akSt ' 

This is the exact solution for the discrete problem with 
St > 0. In the continuum limit, we can take the leading 
order as St — > 0. (| A5|) then reduces to the form The 
recursion for u becomes 

-m 2 u j+1 + Xj + au 2 +l + 0(St), (A7) 



St 
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which is the discrete version of (0. 

Let us now show the connection with the MSR path 
integral discussed in section HlBl We discretize the action 
(|13p with time step St using the ltd prescription. Keeping 
Uj fixed, the path integral formula (|12[) for the generating 
function ([3]) gives us the generating function for Uj+i as 



E{e x ^\u ) 



du 



2tt 



(A8) 



-+m 2 {iij —thj) 5t+u^j rl aiij dt+Xuj-^i 



To leading order for St —> and substituting A — > u this 
becomes identical to the generating function (|A3|) . Note 
that while the first-passage prescription used to obtain 
(|A3|) assumed Uj+i > 0, in (|A8[) we formally allow the 
velocity Uj+i to take any value between — oo and oo. Sur- 
prisingly, this yields the same result to leading order in 
St. It would be interesting to understand how a more 
rigorous MSR approach could be developed directly on 
the discrete version for finite St using first passage times. 



The integrals over Uj+i and Uj+i can be performed ex- 
plicitly, and yield (taking into account St > 0, a > 0, 
w > 0, and iij > 0) 

S(e A " J+1 \uj) = exp [(A - m 2 St + a\ 2 5t) ii, + Xm 2 Stw :i ] . 

(A9) 



Analogously, one can derive a discretized path integral 
for the position variable u for the one-sided Brownian 
potential discussed in section IV Bl 



Appendix B: Functional determinants and 1-loop diagrams 



Here we compute tr(i?E)™ where R is given in (|84p and £ in (|85p . For simplicity we set a = m = 1. Let us recall 
that in Ito discretization 9(0) = 0. First, note that 

(i? T £ T ) tl , t2 = J u(t')sgn(t' - t 3 ) [o(t' - t^'-^ - 9(t 2 - t^*'-^ 

Applying this to tr(iffi)" = tr(i? T £ T ) I \ one gets 



tr(i?E)' 



The convention is that t n +i = t\ and i^+i = tj. Now, we conjecture that for any i , x ...i' rn ^ 

n 

JT sgn(^. - tj ) \6{t' j+1 tj )e-^-^ - 6(t j+1 - t^^A = (-1)" +1 . 



(Bl) 



We were unable to find an analytic proof, but verified this conjucture for n < 5. Assuming it for any n, one obtains 
as claimed 

tr(i?£)™ = - -yfi(t) 
For the one-sided Brownian correlator (|74p we find the self-energy analogous to (|85p as 

E tllta = -25(«! - fe) ^ - h) - 2u(t 2 )9(h - t a ). 

This implies 

(i? T S T ) tl , t2 = -2 / u(i') fff(t a - t')6{t' - ti)e-V-^ + 9(t' - t 2 )9(t 2 - t^e'^-^ 



(B2) 



One then finds tr(iJS) n = for n > I hence a unit functional determinant as claimed in the text. 

This can be generalized to the <i-dimensional interface. We need to compute the functional determinant det(l + i?E) 
with 



Rxiti,x 2 t 2 ~ ^t 1 t 2 (dt 2 S Xl x 2 +32:1x2) 

£ Xl t u x 2 t 2 = S XlX2 [aS tlt2 tixt'Sgnih - t') + au xt2 sgn(t 2 - h)]. 



(B3) 
(B4) 
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We conjecture that this yields 



lndet(l + RE) = trln ^5 XX , - ag xx > J^x'tj = trln 5 XX , - ag xx , J g x , y J^X yt 



(B5) 



For the last equality, we used J t u x t = g xx < J t \ X 't- For a uniform source one recovers the expression in the text of 
section IV Al 

I 



Appendix C: 1-loop functional RG at finite velocity 

In [62|, the 1-loop functional RG equations for a d- 
dimensional elastic interface at non-zero driving velocity 
v > were derived in the Wilson RG scheme. These 
equations have resisted analytical (or numerical) solution 
since then. Here, instead of using Wilson RG with a 
hard cutoff in momentum space, we regularize our model 
by a parabolic well with curvature to 2 . We point out 
that the stationary ABBM disorder correlator (|64[) . (|75l) 
yields a simple solution of the corresponding functional 

I 



RG equations. This also provides an independent check 
of the non-renormalization property for ABBM disorder 
discussed in section IIVBI using a different method. 

For a d-dimensional interface driven by a parabolic well 



of curvature m centered at 



vt, one can derive the 



functional RG flow equation by computing — md m T and 
reexpressing it as a function of T. This is done order by 
order in A, which in this Appendix denotes the renor- 
malized second cumulant of the disorder (the local part 
of the term uu in Y). The resulting functional RG flow 
of A at finite driving velocity v is [63| 



md m A(u) = (e - 2C)A(u) + (uA'(u) 

r°° r°° e -(si+s 2 ) 
+ / dsi / ds2 

Jo Jo s l + s 2 



|A"(m) A(w(s 2 - si)) - A(u + v(s 2 - si)) 
-A'(u + vsi)A'(u-vs 2 ) + A'(u(si + s 2 )) A'(u - vs x ) - A'(u + vs 2 ) } 



(CI) 



Here e = 4 — d, the rescaled correlator is defined via 
A(«) - A d m^A(um() withA" 1 = e J ^(l+k 2 )~ 2 , 
and v = r) m v/m 2 ~ ( ' flows as 

-md m kiv = z-( = 2-(- [ er s A" (sv).(C2) 

Js>0 

The flow of v arises because the friction is corrected by 
disorder. In general, this leads to a non-trivial dynamical 
exponent z defined by the relation above. For v — > one 
recovers the flow at the depinning threshold obtained in 
p2^ . These equations are sufficient (64| for an expansion 
in e with A = O(e). 

Plugging in the correlator for ABBM-type disorder, 
A(u) = A(0) - cr|u|, and C = e into JCl]), one finds 

- md m A{u) = -a 2 (C3) 
-md m v =z-( = 2-e. (C4) 



We see that the dynamical exponent z for ABBM-type 
disorder takes the value z = 2 in any dimension d. The 
ABBM form of the disorder is preserved with —md m a = 
and only A(0) flowing as — md m A(0) = —a 2 . This is 
consistent (for d = 0) with equation (j76|) . In addition, 
as discussed in section IVl and appendix [Bl 2- and higher- 
loop corrections vanish in any d for monotonous motion 
in ABBM-type disorder. More precisely, r — S is the 
logarithm of a functional determinant computed in sec- 
tion El This shows that for ABBM-type disorder, (|CT|) 
is exact to all orders in e = 4 — d. 

We note that for ABBM disorder the correlator re- 
mains non-analytic for any v [65| . This is, presumably, a 
peculiarity of ABBM disorder. For short-ranged disorder 
this may only hold until some scale, the non-analyticity 
being rounded at larger scales (small to). However fur- 
ther studies are needed to clarify the validity of this hy- 
pothesis. 
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